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ABSTRACT 


The  signal  detection  characteristics  of  a  two-input  channel  receiver 
are  studied.  The  signals  received  on  both  channels  are  random  Gaussian, 
narrowband  processes  identical  except  for  a  known  channel-to-channel 
time  delay.  For  narrowband  representation,  the  time  delay  is  treated  as 
a  differential  phase  modulation  of  the  two  signals.  The  noise  on  each  of 
the  two  receiver  channels  is  also  random  Gaussian  processes  which  are 
jointly  wide  sense  stationary  and  correlated.  The  receiver  detects  the 
presence  of  a  signal  by  estimating  the  average  phase  difference  between 
samples  of  the  observed  signals.  The  average  phase  is  an  input  to  a 
matched  filter  detector  where  the  filter  is  matched  to  the  known  phase 
between  the  input  signals.  The  receiver  is  intended  to  detect  at  signal- 
to-noise  ratios  less  than  0  dB  and  over  a  wide  range  of  noise  correlation 
conditions. 

The  theoretical  statistics  required  to  predict  the  detection  per¬ 
formance  are  developed.  The  receiver  operating  characteristic  curves 
are  generated  from  the  likelihood  functions  at  the  matched  filter  output. 
An  example  of  a  possible  signal  phase  function  is  considered  and  the 
receiver  performance  predicted. 

The  performance  of  the  signal  phase  matched  filter  receiver  is 
compared  with  a  cross  correlation  receiver  which  includes  a  noise 
decorrelator  and  signal  phase  compensator.  The  cross  correlation  receiver 
is  more  sensitive  than  the  signal  phase  matched  filter  receiver  if  the 
noise  correlation  is  known.  However,  if  the  noise  correlation  is  unknown 
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and  Che  detection  is  made  using  fixed  signal  processing  and  detection 
threshold,  the  signal  phase  matched  filter  receiver  gives  superior 
performance. 
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CHAPTER  I 


INTRODUCTION 

1.1  Source  of  Detection  Problem 

Radiation  from  a  moving  source  produces  a  varying  time  delay  between 
signals  received  by  two  sensors  which  are  separated  by  a  baseline.  If 
the  geometry  of  the  moving  source  and  the  sensors  is  known,  the  varying 
time  delay  can  be  predicted  prior  to  detection  of  the  signals.  A 
receiver  which  uses  this  known,  varying  time  delay  as  a  basis  to  detect 
the  presence  of  a  signal  source  moving  on  a  known  course  through  a  field 
of  stationary  interfering  noise  sources  has  been  investigated  and  its 
detection  characteristics  have  been  defined. 

Figure  1-1  shows  a  signal  source  which  is  passing  with  a  known  motion 
through  a  field  of  stationary  noise  sources.  The  moving  source  radiates 
a  Gaussian  signal  which  is  received  by  two  physically  separated,  narrow- 
band  sensors.  The  time  delay  between  the  Guassian,  narrowband  signals 
at  the  sensor  outputs  is  determined  by  the  source  motion  and  the  relative 
geometry  of  the  source  and  sensors.  The  time  delay  is  treated  as  a  phase 
delay  in  expressing  the  received,  narrowband  signals. 

The  presence  of  stationary  interference  noise  sources  in  an  isotropic 
Gaussian  background  noise  field  results  in  an  anistropic  noise  field. 

This  anistrophic  noise  combines  with  the  electrical  noise  of  the  sensors 
such  that  the  noise  at  the  filtered  sensor  outputs  is  correlated.  The 
sensor  output  noise  is  stationary  and  Gaussian  with  unknown  power  level. 

This  research  studied  the  detection  problem  after  the  reception  of 
the  signals  and  noise  by  the  two  sensors.  The  problem  investigated  is 
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the  signal  detection  characteristics  of  a  receiver  with  two  input  channels. 
The  signals  in  both  channels  are  narrowband,  random,  Gaussian  processes 
which  are  identical  except  for  a  known,  time  varying  phase  delay.  The 
input  noise  is  correlated,  Gaussian  noise  of  unknown  power  level.  The 
receiver  detection  process  estimates  the  average  phase  difference  between 
the  observed  signals  in  the  two  receiver  input  channels.  The  estimated 
average  phase  is  an  input  to  a  detector  which  is  an  approximation  of  a 
matched  filter  detector  for  the  known  phase  delay  between  the  desired 
input  signals. 

1.2  Assumptions  and  Limitations 

The  signal  detection  system  is  intended  to  function  at  low  signal- 
to-noise  ratios  and  over  a  variety  of  noise  correlation  values.  The 
signal  processing  is  to  be  implemented  using  digital  techniques.  These 
factors  and  the  nature  of  the  original  detection  problem  result  in  a 
number  of  assumptions.  Any  limitations  imposed  by  these  assumptions 
must  be  considered  when  interpreting  the  results. 

The  received  signal  is  a  stationary  narrowband,  zero-mean  Gaussian 
stochastic  process.  The  signal  in  one  input  is  a  phase-shifted  version 
of  the  signal  in  the  other  input.  The  signal  power  level  is  equal  in 
each  input.  The  time  varying  phase  delay  between  the  input  signals  is 
known.  The  bandwidth  of  the  input  signals  is  much  greater  than  the  band¬ 
width  associated  with  the  time  varying  phase  delay. 

The  noise  in  the  two  inputs  is  jointly  stationary,  correlated,  zero- 
mean,  narrowband  Gaussian  noise.  The  power  level  of  the  noise  is  unknown, 
but  it  is  the  same  in  each  input.  The  detection  process  is  investigated 
for  both  known  and  unknown  noise  correlation  cases. 


The  system  is  intended  to  detect  at  signal-to-noise  ratios  less 
than  0  dB.  The  signal-processing  interval  is  selected  to  be  long  enough 
to  achieve  this  goal. 

The  inputs  are  sampled  and  digitized.  It  is  assumed  that  the 
sampling  period  is  adjusted  to  the  input  bandwidth  such  that  sequential 
samples  are  statistically  independent. 


1. 3  Background 

Detection  of  a  Gaussian  signal  in  Gaussian  noise  has  been  widely 
investigated  in  the  past  [1],  [2],  [3].  The  likelihood  detector  is 
commonly  used  to  detect  signals  which  are  not  completely  known.  If  the 
input  noise  is  correlated,  the  likelihood  detector  can  be  implemented 
using  a  noise  de-correlator  or  spatial  prewhitener  [4]*  These  investiga¬ 
tions  assumed  the  noise  power  level  and  noise  correlation  were  known. 

If  the  power  level  is  unknown,  a  detector  can  be  implemented  consisting 
of  two  channels  [5],  [6].  One  channel  has  the  form  of  a  likelihood  ratio 
detector  and  the  other  channel  is  a  minimum  variance  power  level  estima¬ 
tor.  With  two  inputs,  this  detector  takes  the  form  of  a  cross  correla¬ 


tion  detector.  Source  motion  causes  degradation  in  the  output  of  the 


correlator,  unless  the  correlator  implements  a  compensation  delay 
modulation,  [7],  [8]. 

Based  on  the  above  background,  a  receiver  which  could  be  used  to 
solve  the  detection  problem  would  consist  of  a  noise  decorrelator,  time 
delay  compensation  modulator  and  cross  correlation  detector.  Such  a 
receiver  requires  that  the  noise  correlation  be  known  in  order  to  obtain 


maximum  sensitivity. 
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An  alternative  to  the  classical  likelihood  detector  is  a  receiver 
which  estimates  the  phase  between  the  signals  at  the  filtered  sensor 
outputs.  The  signal  is  then  detected  using  a  matched  filter  detector 
[9].  The  filter  is  matched  to  the  known  phase  delay.  This  alternate 
receiver  is  analyzed  and  its  detection  characteristics  determined. 

1.4  Organization 

Chapter  II  contains  a  detailed  description  of  the  detection  problem. 
The  generalized  received  signal  and  noise  model  is  explained.  This  is 
followed  by  a  description  of  the  methods  of  estimating  phase  difference 
and  a  description  of  the  matched  filter  detector. 

The  theoretical  statistics  required  to  determine  the  detection 
system  performance  are  developed  in  Chapter  III.  These  statistics  are 
a  function  of  the  input  signal-to-noise  ratio  and  the  correlation 
characteristics  of  the  noise  in  the  two  channels  of  the  receiver.  Major 
emphasis  is  placed  on  the  statistics  of  the  average  phase  estimate.  The 
material  in  this  chapter  forms  the  theoretical  basis  for  understanding 
the  signal  detection  problem. 

In  Chapter  IV,  the  receiver  detection  characteristics  are  determined 
for  a  linear  phase  delay  function.  The  matched  filter  configuration  is 
described,  the  effect  of  phase  averaging  examined  and  the  theoretical 
detection  curves  generated.  The  above  assumes  the  noise  correlation  is 
known.  The  chapter  ends  by  investigating  detection  performance  for  the 
unknown  noise  correlation  case. 

A  performance  comparison  with  the  cross  correlation  detector  is 


made.  This  performance  comparison  is  made  for  both  known  noise 
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correlation  and  unknown  noise  correlation  cases.  The  cross  correlation 
detector  with  noise  de-correlator  and  phase  shift  compensation  is 
described  in  Appendix  B. 

Conclusions  and  recommendations  are  contained  in  Chapter  V. 


CHAPTER  II 


DETECTION  PROBLEM 

Detecting  a  dual  channel  differential  phase  modulated  signal  in 
the  presence  of  correlated  noise  is  a  problem  in  selecting  one  of  two 
possible  hypotheses,  the  desired  signal  is  present  or  the  desired  signal 
is  not  present.  The  receiver  being  studied  uses  the  Neyman-Pearson 
criterion  [10]  to  test  the  hypotheses.  This  criterian  maximizes  the 
probability  of  detection  for  a  given  probability  of  false  alarm.  The 
observable  signals  at  the  receiver  inputs  may  be  of  two  forms: 


r-^t) 


+  nL(t) 


r2(t) 


+  n2(t) 


where  r(t)  =  observable  signal 
s(t)  =  desired  signal 
n(t)  =  noise 
x(t)  =  time  delay. 

The  receiver  is  designed  to  operate  on  the  observable  signals  and  select 
either  the  null  hypothesis,  H^,  or  the  alternative  hypothesis,  H^. 
Symbolically  these  hypotheses  are: 


H  : 
o 


r 1 ( t)  =  n1(t) 
r2(c)  =  n2(t) 


,  0  £  t  £  T 
,  0  £  t  £  T 


(2-1) 
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[  r^t)  =  s(t)  +  nL(t)  ,  0  <  t  *  T 

H  j  / 

1  |  r2(t)  =  s^t-Tg(t)j  +  n2(t)  ,  0  5  t  i  T 

where  T  =  period  of  s(t). 

Figure  2-1  shows  a  block  diagram  of  the  receiver  being  studied. 

The  receiver  consists  of  four  basic  blocks.  The  observed  signals  are 
sampled  and  digitized.  An  estimate  is  made  of  the  average  difference 
between  phases  of  the  two  observable  sequences.  The  resulting  average 
phase  sequence  is  an  input  to  a  matched  filter.  The  transfer  function 
of  the  filter  is  matched  to  the  differential  phase  shift  between  the 
desired  signals  received  in  the  two  input  channels.  If  the  filter  output 
exceeds  a  threshold,  the  hypothesis,  signal  present,  is  selected.  If 
the  threshold  is  not  exceeded,  the  Hq  hypothesis,  no  signal  present,  is 
selected. 

2.1  Signal  and  Noise  Model 

The  signal-noise  model  for  the  generalized  detection  problem  is 
shown  in  Figure  2-2.  The  noise  sources  n^Ct) ,  nv(t)  and  n^(t)  are  always 
present.  The  signal  s(t),  which  may  or  may  not  be  present,  is  a  random 
narrowband  Gaussian  process  with  zero  mean.  A  random  process  is  narrow- 
band  if  its  spectral  density  is  zero  except  for  a  narrow  region  around 
a  high  carrier  frequency  [11].  The  received  signal  can  be  expressed  as: 

s(t)  =  s  (t)  cos  ui  t  -  s  (t)  sin  w  t 
c  c  s  c 


(2-3) 
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where  sc(t)  *  the  in-phase  component  of  s(t) 

sg(t)  =»  the  quadrature  component  of  s(t) 

■  the  center  frequency  of  the  narrowband  spectral 
density  of  s(t) . 

It  is  assumed  that  the  maximum  rate  with  which  the  time  delay,  t  (t) , 

s 

varies  is  very  small  when  compared  to  the  inverse  of  the  signal  band¬ 
width.  The  signal  bandwidth  is  defined  by  the  receiver  passband.  The 
in-phase  and  quadrature  components  of  the  time-delayed  signal  can  there¬ 
fore  be  considered  to  be: 


sc[t-i:s(t)]  -  Sc(t) 

Ss[t_Ts(t)]  -ss(t)* 

The  time-delayed  signal  is  approximated  as: 

s  [t“ Ts(t)j  -  sc(t)  cos  [a)ct-ucTs(t)]  _  sg(t)  sin  [V-Vs<‘>]  ’ 

If  the  u  x  (t)  is  considered  to  be  a  variable  phase, 
c  s 

*s(t)  =  ucTs(t)’ 

the  time  delayed  signal  may  be  expressed  as: 

s£t-Tg(t)j  =  Sc(t)  COS  [uct-!t,g(t)]  ~  Sg(t)  sin  ’ 

The  signals  in  the  two  channels  of  the  receiver  are  expressed  as: 

s..(t)  =  s  (t)  cos  a)  t  -  s  (t)  sin  w  t 
i  c  c  s  c 

s2(t)  =  sc(t)  cos  [»ct-*a(t)j  -  ss(t)  sin  ^0)ct-^s(t)j  .  j 


(2-4) 


(2-5) 
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The  total  received  noise  is  the  sum  of  the  contributions  from  the 
receiver  electrical  noise,  the  isotropic  noise  field  and  all  stationary 
interference  sources.  The  received  noise  is  represented  by  three 
independent  noise  sources  in  the  signal-noise  model.  The  noise  in  each 
channel  can  be  considered  to  be  the  sum  of  a  common  or  fully  correlated 
component  and  an  independent  component.  The  independent  components  are 
represented  by  n^(t)  and  n^Ct).  The  common  component  is  n^(t).  The 
fixed  time  delay  rn  associated  with  n^(t)  is  a  composite  of  the  time 
delays  of  all  the  interference  sources. 

As  these  noise  sources  are  independent  random  narrowband  processes, 
they  can  be  expressed: 


n  (t)  =  n  (t)  cos  to  t  -  n  (t)  sin  10  t 
u  uc  c  us  c 

n  (t)  =  n  (t)  cos  ui  t  -  n  (t)  sin  to  t 
v  vc  c  vs  c 

n  (t)  =  n  (t)  cos  to  t  -  n_  (t)  sin  to  t 
I  Ic  c  Is  c 

n^t— rn)  =  nIc(t)  cos  (w  t-^)  -  nIs(t)  sin  (to  t-^) 


)  (2-6) 


I 


where 


<o 


T 

c  n 


Using  the  above  narrowband  representation,  the  noise  in  the  two  receiver 
channels  is : 


nl(t)  =l 

n  (t) 

.  uc 

+  nIc(t>- 

cos  to  t  -  [n  (t)  +  nT  (t)  1  sin  to  t 
|  c  [us  Is  J  c 

n2(t)  = 

n  (t) 

.  vc 

+  nlc(t) 

cos  +  n^g(t)  sin  j  cos  toct 

>  (2-7) 

-  [nvs(C)  “  n^c(t)  S*n  ^1  +  nIs^  COS  ^1  sin  0Jct  ’  j 
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The  in-phase  and  quadrature  components  of  n^(t)  and  n^(t)  are: 


nlc(t)  =  nuc(t)  +  nIc(t) 


nls(t)  =  nus(t)  +  nIs(t) 


n2c(t)  =  nvc(t)  +  nIc(t)  COS  +  nIs(t)  sin  *1 
n2s(t)  *  nvs(t)  "  nIc(t)  Sin  +  nls(t)  cos  *i 


>  (2-8) 


The  mean,  variance  and  covariance  of  the  in-phase  and  quadrature 
components  of  the  noise  are  important  parameters  in  determining  the 
detection  statistics.  It  is  assumed  that  the  receiver  is  designed  such 
that  the  power  level  of  the  noise  in  each  channel  is  equal. 


2  2 
a  =  a 
u  v 


(2-9) 


Since  the  noise  sources  are  all  assumed  to  be  zero  mean  processes,  the 
mean  values  of  the  noise  components  are: 


E{nlc(t)}  =  E{ni3(t)}  =  E{n2c(t) }  *  E{n2s(t)  }  =  0  ; 


(2-10) 


the  variances  are: 


E(nlc2  (t)>  -  au2  +  ax2  =  a,2 


2  .  .  ,  2  ,  2  2 
E{n  (t) }  =  a  +  oT  =  o„ 
Is  u  I  N 


E{n2c2  (t)  }  =  c^2  +  Oj2  (cos2  $x  +  sin2  $  ) 
E{nXg2  (t)}  =  c^2  +  a^.2  (sin2  +  cos2  $  ) 


=  a.. 


=  a  . 


>  (2-11) 
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the  covariances  become 
Mlc,ls  “  y2c,2s  =  0 

2  2 

ulc,2c  *  uls,2s  "  PN°N  =  °N  KNC0S4lN 

2  2 

“ylc,2s  “  uls,2c  “  XN0N  =  CN 

Thus,  the  noise  in  the  two  channels  of  the  receiver  is  correlated. 

The  observable  signals  are  the  sum  of  equations  (2-5)  and  (2-7). 

Using  the  noise  components  defined  by  equations  (2-8),  the  observable 
signals  become: 

rl(t)  =  [sc(t)+nlc(t) ]  cos  (i>ct  -  £ss(t)+nls(t)]  sin  u)ct 

r2(t)  =  [sc(t)  cos  4>s(t)+sg(t)  sin  <{>g(t)+n2c(t) j  cos  u^t  >  (2-15) 

-  [ -s  (t)  sin  <j>  (t)+s  (t)  cos  <f>  (t)+n~  (t)|  sin  u  t 
L  c  s  s  s  <*s  j  c 


where  r(t)  is  the  Hilbert  Transform  of  r(t).  For  a  narrowband  signal 
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r(t)  the  analytic  signal  is: 


rp(t)  =[rc(t)  +  j  rs(t)]  eJ  C 

The  analytic  signal  is  sampled  to  obtain  the  complex  low  pass  signal 
r^(t)  +  j  rg(t).  Thus,  the  sampled  signal  is: 

rp(n)  =  rc(n)  +  j  rs(n) 

The  sampling  is  at  a  rate  which  is  equal  to  or  greater  than  the  inverse 
of  the  bandwidth  of  the  low  pass  signal  [13]. 

After  quadrature  sampling,  the  observable  signals  can  be  expressed 


as : 


rpl(n) 


rP2(n) 


[sc(n)  +  nlc(n)]  +  j[Ss(n)  +  nls(n)] 

|s  (n)  cos  4)  (n)  +  s  (n)  sin  <j>  (n)  +  n  (n) 

L  c  s  s  s  Zc  J 

+  j [-sc(n)  sin  $g(n)  +  sg(n)  cos  4g(n)  +  n2g(n)j 


}  <2 


The  complex  sequences  rp^(n)  and  r^Cn)  ma^  exPresse^  in  terms  of 
their  amplitude  and  phase, 


rpl(n)  =  A^n)  e 


(2 


where 


Ai(n)  =  {  [sc(n)+nlc(n)j 2  +  ^sg(n)+nlg(n)j 2  } 


1/2 


i  (n)  =  tan 


-1 


s  (n)  +  n.  (n) 
s _ is 

sc(n)  +  r>lc(n) 


and 

j^2(n) 

rp2(n)  =  A2(n)  e 


( 


-16) 


-17) 


-18) 


where 


A9(n)  =  { fs  (n)  cos  $  (n)  +  s  (n)  sin  $  (n)  +  n9  (n)l^ 

—  L  C  s  s  s  —  c  J 

r  i  2  ^ 

+  -s  (n)  sin  a  (n)  +  s  (n)  cos  $  (n)  +  n„  (n)  } 

L  c  s  c  s  2s  J 


$9(n)  =  can 


,  -s  (n)  sin  <4  (n)  +  s  (n)  cos  i>  (n)  +  n.  (n) 
-1  c _ _s _ c _ s _ 2s 

s  (n)  cos  $  (n)  +  s  (n)  sin  A  (n)  +  n„  (n) 
*•  c  s  s  s  2c 


2. 3  Average  Phase  Estimate 

It  has  been  noted  that  the  bandwidth  of  the  received  signal  is 

assumed  to  be  much  greater  than  the  bandwidth  of  the  time  delay  function. 

Therefore,  processing  of  the  estimated  phase  difference  between  r^(n) 

and  r  „(n)  mav  be  simplified  by  averaging.  Two  methods  that  can  be  used 

p2 

to  make  an  average  phase  estimate  are  shown  in  Figure  2-4.  These  two 
methods  are  identified  as  angle  average  and  vector  average. 

If  the  complex  conjugate  of  rp2(n)  i-s  multiplied  by  rp^(n)  the 
product  will  be: 


e(n)  =  rpl(n)  r?2  (n) 


A^(n)  A2(n)  e 


j  ^1(n)-<J>2(n)j 


(2-19) 


Thus,  the  phase  of  e(n)  is  the  phase  difference  of  sequences  r^(n)  and 
rp2(n).  Assuming  that  no  noise  is  present,  the  signal  phase  in  each  of 
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the  receiver  channels  is: 


4^(n)  =  can 


s  (n)  I 
s 

s  (n) 
r 


4_,(n)  =  tan 


s  (n)  cos  4  (n)  -  s  (n)  sin  4  (n) 
s _ _s _ c _ s 

s  (n)  cos  4  (n)  +  s  (n)  sin  4  (n) 
c  ns  s 


(2-20) 


s  (n) 
n 

s  (n) 
c 


r  sin  4  (n)  ^ 

I  - § - 

[  cos  4>s(n) 


s  (n)1f  sin  4  (n) ' 


T  s  (nnr 


cos  4  (r 
s 


Since 


/  _  tan  a  -  tan  8 

tan  (a-S)  =  - -  , 

1  +  tan  a  tan  a 

the  phase  in  channel  2  is 


|>2(n)  =  41(n)  -  4g(n) 


(2-21) 


The  phase  difference  is  then 


j>,(n)  -  (n)  =  4  (n) 

i  1  s 


(2-22) 


The  phase  of  e(n)  is  the  difference  between  the  phase  of  the  signals 
in  the  two  receiver  channels.  When  noise  is  added  to  the  observable 
signal,  the  assumption  that  the  phase  of  e(n)  is  the  desired  signal 
phase  difference  is  no  longer  certain. 

To  estimate  the  average  phase  of  the  complex  sequence  e(n),  the 
sequence  is  first  divided  into  its  real  and  imaginary  components.  If 


,  Vn> 

e(n)  =  [e(n)  [  e-1 
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2. 4  Matched  Filter  Detector 

The  matched  filter  transfer  function  is  defined  to  be: 

H(ju>)  -  S*(jw)  e~2“T  (2-27) 

where  S(joj)  =  the  Fourier  transform  of  the  desired  signal 
T  =  processing  interval. 

If  the  signal  has  a  period  T  and  is  sampled  and  expanded  in  a  Discrete 
Fourier  Series  [14),  the  digital  equivalent  to  the  matched  filter 


22 


23 


Expressing  equation  (2-31)  using  the  sine  and  cosine  Fourier  components 


V1 


»T  *  %<»>  V0)+  5  [Vk)+JVH  [*s(W-ib,(M] 


XT  =  a  (0)  a  (0)  + 

L  $  S 


V1 


1  V1 
I  Z  [ 

k=l  <■ 


a  (k)  a  (k)  +  b  (k)  b  (k)l 
$  s  $>  s  J 


1  Z  k<» 

F  k=l  L 


b  (k)  -  a  (k)  b  (k)l 

<$>  <t>  s  J 


For  a  real  sequence,  the  cosine  terms  are  even  functions  and  the  sine 
terms  are  odd  functions 
a(k)  =  a(n-k) 
b(k)  =  -b(n-k) 

Therefore,  the  output  of  the  matched  filter  can  be  expressed  as 


-1 


XT  =  a  (0)  a  (0) 

L  <p  S 


+  Y  [a  (k)  a  (k)  +  b  (k)  b  (k)l 
NF  k=l  L  $  s  $  s  J 


nf  nf 

+  ax(~?)  a  >  N  even 

(J)  l  S  l 


If  the  a  (0)  and  a  (^r)  terms  equal  zero,  the  output  of  the  matched 
s  si 


filter  is 


XT  = 


-1 


T  Z  [ 

F  k=l  L 


_  .  .  . »  N  e\ 

S  s  J 


a  (k)  a  (k)  +  bA(k)  b  (k)|  ,  N  even 


(2-32) 


V1 


(  4  [%(k)  as(k)  +  Vk)  bs(k)]  ’ 


N  odd 
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The  common  multiplier  of  two  has  been  dropped  since  it  applied  to  all 

N 

terms  other  than  a(0)  and  aOj) •  The  two  versions,  complex  coefficients 
and  sine-cosine  coefficients ,  of  a  discrete  time  matched  filter  are 
shown  in  Figure  2-5. 

The  output  of  the  matched  filter  is  the  input  to  the  threshold 
detector.  The  threshold  is  selected  for  a  given  probability  of  false 
alarm.  If  the  output  exceeds  the  threshold,  hypothesis  is  selected 
and  the  signal  is  considered  to  be  present. 
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CHAPTER  III 

THEORETICAL  DEVELOPMENT 


In  order  to  establish  the  receiver  performance,  it  is  necessary  to 
know  the  statistics  at  the  threshold  detector.  In  this  chapter,  the 
statistical  properties  of  the  random  processes  are  traced  through  the 
signal  processing  system  shown  in  Figure  2-1.  The  statistical  properties 
of  the  observable  sequences  are  defined  from  the  signal  and  noise  model. 

The  mean  and  variance  of  the  average  phase  estimate  is  calculated  in 
terms  of  the  signal-to-noise  ratio,  signal  phase,  and  of  the  complex 
correlation  coefficient  of  the  noise  present  in  the  observable  sequences. 

The  statistics  of  the  Fourier  coefficient  are  determined  for  a  non-stationary 
average  phase  estimate.  From  the  statistics  of  the  Fourier  coefficient,  the 
statistics  at  the  threshold  detector  are  defined  and  the  probability  of 
false  alarm  and  probability  of  detection  determined. 


3.1  Input  Signal  and  Noise  Characteristics 

The  statistical  properties  of  the  observable  complex  sequence 

r  , (n)  and  r  »(n)  given  in  equations  (2-16)  can  be  defined  in  terms  of  the 
pi  p/ 

signal  and  noise  model  described  in  Chapter  II. 

The  signal  is  a  wide-sense  stationary,  narrowband,  zero-mean, 

2 

Gaussian  stochastic  process  having  a  power  of  a  .  The  mean  and  variance 

s 

of  the  in-phase  and  quadrature  components  of  the  signal  are: 


E{s  (t) }  =  E{s  (t) }  ■  0 
c  s 

E{s  2(t)}  =  g{s  2(t)}  =  o 

c  S3 


(3-1) 


2 


(3-2) 
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E{s  (t)  s  (t)}  =  0  (3-3) 

c  s 

The  received  noise  is  correlated,  wide-sense  stationary,  zero-mean, 

2 

narrowband  Gaussian  noise  with  power  a..  .  The  mean  and  covariance  of  the 

N 

in-phase  and  quadrature  components  of  the  noise  are  given  by  equations 
(2-10)  through  (2-14). 

The  observable  complex  sequences  can  be  defined  in  terms  of  the 
in-phase  and  quadrature  components . 


From 


form 


rpl(n)  =  rlc(n)  +  j  rls(n) 


rP2(n)  =  r2c(n)  +  j  r2s(n) 


Equation  (2-15) : 


rlc<n)  =*  s£(n)  +  Nlc(n) 


rls(n)  =  Ss(n)  +  Nls(n) 


(3-4) 


(3-5) 

r9  (n)  ■  s  (n)  cos  <p  (n)  +  s  (n)  sin  <j>  (n)  +  n  (n) 

ZC  C  S  S  S  tC 

r2s(n)  =  sc(n)  sin  4>s (n)  +  Sg(n)  cos  <pg( n)  +  N2g(n) 

The  in-phase  and  quadrature  components  of  the  observable  sequences 
a  vector. 


» 

I  =  [rlc(n),rls(n) ,r2c(n) ,r2s(n) ]  (3-' 

with  a  four-dimensional  Gaussian  probability  density  function 

p(r  ,r,  ,r„  ,r.  ).  The  covariance  matrix  of  the  four-dimensional  proba- 
lc  Is  2c  2s 

bility  density  function  can  be  shown  to  be: 
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R  =  a 
r  N 


1  +  h 


1  +  h 


p,.+h  cos  6  A.,+h  sin  4> 

.M  S  N  S 


-X„-h  sin  it  p„+h  cos  d> 

N  s  N  s 


p+h  cos  4>  -A..-h  sin 

M  s  N 


s' 


A  +h  sin  <j)  P„+h  cos  <p 

N  s  N  s 


1  +  h 


1  +  h 


(3-7) 


where  h  =  — ~  ,  signal-to-noise  ratio. 
°N 


(3-8) 


The  covariance  matrix  is  positive  definite.  The  four-dimensional  Gaussian 
probability  density  function  is: 


p(r_)  =  - = — 

(2  )  |R 


,  1  ’  -1  , 
-q  exP  (~~2—  Rr  i) 


(3-9) 


where 


R  —  determinant  of  matrix  R 
1  r 1  r 

R  ^  —  inverse  of  matrix  R  . 
r  r 


The  correlation  coefficient  of  the  sequence  is 


K  e 
o 


j*o  A 


E{r1(n)r;,  (n)} 


^E{r1(n)r1*(n) }  E  (r2(n)r2  (n)}j 


(3-10) 


Substituting  the  appropriate  terms  from  the  covariance  matrix,  the 


correlation  coefficient  becomes 
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K 


J9o  (h  cos4>s  +  +  j  (H  sin£s  +  AN> 


1  +  h 


(3-11) 


The  magnitude  and  phase  of  the  correlation  coefficient  are: 


K 


2  2  2 

~  h  +  p  +  A  +  2ph  cos  P  +2  A  h  sin  <p 

_  N  N  vl  S  S 


tan 


-i 


(1  +  h) 

h  sin  tp  +  A„ 
s _ N 

h  cos  0  +  pM  : 

s  N  j 


(3-12) 


(3-13) 


3.2  Statistical  Properties  of  the  Average  Phase  Estimate 

The  processes  for  obtaining  the  average  phase  estimates  have  been 

described  in  Chapter  II.  The  average  phase  statistics  will  be  shown  to 

be  a  function  of  the  complex  correlation  coefficient  between  the  two 

observable  sequences  r  . (n)  and  r  „(n). 

pi  p2 

The  phase  difference  between  rpl(n)  and  r^2(n)  has  a  probability 
density  function  given  by  [15] 


1-K 


p(4>)  =  '5°  f — +  — 2  3/2  +  Sin_1  B)  } 

*  1-8  (1-3  )' 


(3-14) 


where 


0  =  K  cos  (4>— 4>  ),  -  -^  <  sin  1  3—  -r 

O  0  4  —  4 


(3-15) 


The  mean  and  variance  of  the  single  point  phase  estimate  are  [16] 


cos  ^(KqcosP0) 

v7T"2  27 

1-K  cos  P 

O  0 


K  sin  p 
o  0 


(3-16) 
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a  2 

<P 


12  ~  E(Ko>  + 


-1 

cos 


(K  cos4>  ) 
o  o 


2_oV>_ 

(1-K  2cos2<i  ) 
o  o 


where 


E(K  ) 
o 


^  I 

nal 


o 


(2n)2 


(3-17) 


(3-18) 


The  mean  and  standard  deviation  of  the  single  phase  estimate  are  plotted 

as  a  function  of  K  and  <j>  in  Figure  3-1. 
o  a 

3.2.1  Angle  Average  Phase  Estimate 

The  angle  average  phase  estimate  was  defined  by  Equation  (2-25). 
— .  2 

If  <J(n)  and  (n)  are  constant  over  the  averaging  interval,  the  mean 
and  variance  of  the  angle  average  phase  estimate  are  [17] 


4> 


aa 


2 

a 

aa 


<Kn) 

y<»> 

N 


(3-19) 

(3-20) 


These  values  are  also  in  Figure  3-1.  The  mean  of  the  angle  average  phase 
estimate  is  independent  of  the  number  of  points  averaged.  Thus  Figure  3.1(a) 
also  represents  the  mean  of  angle  average  phase  estimate  plotted  as  a 
function  of  Kq  and  4>q.  Figure  3.1(b)  can  be  considered  as  a  normalized 
plot  of  the  standard  deviation  of  the  angle  average  phase  estimate  by 

labelling  the  absissa  NO  /it. 

SL  dicL 

It  is  assumed  that  the  number  of  points  averaged  is  large  enough 
so  that  the  central  limit  theorem  may  be  invoked.  This  means  that  the 
angle  average  phase  estimate  can  be  assumed  to  have  a  Gaussian  distribution. 
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3.2.2  Vector  Average  Phase  Estimate 

The  vector  average  phase  estimate  was  described  in  Chapter  II 

and  defined  by  equation  (2-26).  Since  an  expression  for  the  probability 

density  function  of  the  vector  average  phase  estimate  is  not  known,  it  is 

necessary  to  derive  such  an  expression.  The  mean,  variance  and  correlation 

coetiicient  are  derived  for  the  real  and  imaginary  components  of  the  complex 

signal  obtained  when  multiplying  complex  sequence  r  , (n)  by  the  complex 

pi 

conjugate  of  the  complex  sequence  r^Cn).  The  real  component  is  considered 
to  be  the  X  vector  and  the  imaginary  component  is  the  Y  vector.  The 
X  and  Y  vector  components  are  then  averaged.  It  is  assumed  that  the  number 
of  points  averaged  is  sufficient  to  invoke  the  central  limit  theorem  and 
consider  the  average  X  and  average  Y  vectors  to  be  Gaussian  distributed. 
These  variables  are  shown  in  Figure  2-4. 

An  expression  is  then  derived  for  the  probability  density  function  of 
the  vector  average  phase  -ssuming  the  average  X  and  average  Y  vectors 
are  correlated  Gaussian  parameters  with  unequal  means  and  variances. 

Since  the  expression  for  the  probability  density  function  is  not  easily 
integrable,  curves  of  the  mean  and  standard  deviation  of  the  vector 
average  phases  are  generated  using  numerical  techniques. 

Expressions  for  the  mean,  variance  and  covariance  of  the  real  and 
imaginary  components  of  the  product  of  a  complex  variable  and  the  complex 
conjugate  of  another  complex  variable  are  derived  in  Appendix  A.  These 
expressions  are  given  in  terms  of  the  components  of  the  covariance  matrix. 
Substituting  the  covariance  matrix  terms  from  equation  (3-7),  the  equations 
for  the  moments  of  the  X  and  Y  vectors  become 
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X  =  2  a  2  K  (1+h)  cos  5 
Me  o 

Y  =  2  a  2  K  (1+h)  sin  4 
No  o 


a  2  =  2 

X 

4 

aN 

(l+h)2(l  + 

x  2 

O 

cos 

2V 

a  2  =  2 
y 

4 

°N 

(l+h)2(l  - 

K2 

o 

cos 

2V 

y  =2 

xy 

4 

2  2 
(1+h)  K 

o 

sin 

2$ 

0 

The  average  X  and  average  Y  vectors  are  calculated  from 

N  -1 
a 

Xa  -  r  I  x(n) 

a 

n=o 


Y 

a 


N  -1 
a 


n=o 


The  moments  of  X  and  Y  are 
a  a 


X  =  X 

a 

Y  =  Y 

a 


(3-21) 

(3-22) 

(3-23) 

(3-24) 

(3-25) 


(3-26) 


(3-27) 


(3-28) 

(3-29) 

(3-30) 

(3-31) 
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U 

a 


(3-32) 


a 

Knowing  the  moments  of  the  real  and  imaginary  vector  components, 
the  vector  average  phase  probability  density  function  can  now  be  derived. 
Assume  the  X  and  Y  vectors  are  jointly  Gaussian  distributed,  the  joint 
probability  density  function  is  given  by: 


P 

xy 


(x,y) 


2to  a  (1-p  ) 

xa  ya  xy" 


2<1-PW  )  L 


«a-V2 

a  2 

xa 


(X-X  ) (Y-Y  ) 

3  3  3  3 

a  o 
xs  va 


(3-33) 


The  X  and  Y  vectors  are  related  to  the  phase  angle,  $  and  amplitude, 
a  a  a 

r  by : 
a  J 


r  =  (X2  +  Y2)h 
a  a 


=  tan 
va 


•1  Ya 


(3-34) 

(3-35) 


a 

The  probability  density  function  of  the  phase  angle,  $va>  can  be 

determined  by  finding  the  joint  probability  density  function  of  r  and 

$  and  then  integrating  with  respect  to  r.  The  joint  probability  densitv 
va 

function  of  r  and  <p  can  be  found  from  p  (x,y) 
a  va  xy 


p(r’Ka} 


thp 


(r  cos  <t> 
xy  va 


r  sin 


$  ) 

va 


(3-36) 


where 
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9r  3r 

3  X  3  Y 


3  tf>  3  4> 

va  va 

3  X  3  Y 
a  a 


(3-37) 


Substituting  Equation  (?-34)  and  (3-35)  in  (3-37),  the  Jacobian  becomes: 


(3-38) 


/  2  2 

X  +  Y 
a  a 


The  joint  probability  density  function  of  r  and  $  then  becomes: 


p”,1^  eiq>  H'-'Wj 

xa  ya  xy 


(3-39) 


where 


,  (r  cos<{>  -X  )  (r  cos<j>  -X  )(r  sm<U  -Y  ) 

/  ,  ,  1  va  a  va  a'  va  a  , 

a(r,$  )  =  -r  - - = - 2p  - - - - - + 

va  2  „  2  xy  o  a 

La  xa  ya 

xa 


(r  sintj)  -Y  )2 
va  a 


(3-40) 


The  probability  density  function  P($va)  can  1)6  found  from 


P(^va)  * 


00 

“7 - j-qJ  r  exp[-a(r,4»va)] 


(3-41) 


2mj  a  (l-p 
xa  ya  xy 


Performing  the  indicated  integration  and  substituting  the  moments  of 


X  and  Y  ,  the  vector  average  phase  probability  densitv  function  becomes 
a  a 
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P^va5 


2ir[l-K  2  cos(4>  — <+>  )] 

o  va  o 


exp  - 


1  +  K 


vfi  K  (1-K  2)  cos  (<fi  -<t>  ) 
a  o _ o_ _ va  o 

2i/tt  [1-K  2  cos(4>  -4>  )]2/*2 
1  r,  v  wa  o 


exp  - 


N  K  2  sin2 (4>  -4>  ) 


1-K  cos  2  (<})  -<J>  ) 

o  va  o 


X  <  1  +  e  r  f 


r  ✓’N  K  (l-K2)1*  cos  (4>  -<J)  )  , 

a  o  va  o 

.  (1+K  2)i5 1 1-K  2  cos  2(4>  -4>  )]  15  . 

o  l  o  va  01 


(3-42) 


where 

=  number  of  points  averaged 

Kq  =  magnitude  of  correlation  coefficient  of  the  receiver  observable 
sequences  r^(n)  and  rp2^n)  complex  sampling 

4>q  =  phase  of  correlation  coefficient 


and 

x 

erf  (x)  =  —  j  exp  (-Z2)  dZ 

/rrJ 

o 

Probability  density  functions  for  the  vector  average  phase  are  shown 
in  Figure  3-2 (a),  for  various  values  of  Kq  when  the  number  of  points 
averaged  is  held  constant,  and  in  Figure  3-2 (b),  for  various  values  of  N 

dk 

when  K  is  constant.  In  both  figures  p(4>  )  is  plotted  as  a  function  of 

O  V3 

(4>  ~4>  )  where  4)  is  the  phase  angle  of  the  correlation  coefficient  defined 

va  0  0 

by  Equation  (3-20).  Since  the  expression  for  the  probability  density 

function  is  not  readily  integrable,  computer  programs  were  developed  to  give 

curves  of  the  mean  and  variance  of  the  vector  average  phase  estimate  as  a 

function  of  4>  a  given  K  and  N  .  Curves  of  the  mean  and  standard 

o  o  a 


deviation  of  the  vector  average  phase  are  shown  in  Figures  (3-3)  to  (3-6). 


FIGURE  3-3.  Mean  and  Standard  Deviation  of  Vector  Average  Phase  Estimate,  K  =  0.025, 


<AI  MEAN  <BJ  STANDARD  DEVIATION 


FIGURE  3-5.  Mean  and  Standard  Deviation  of  Vector  Average  Phase  Estimate,  K..  =  0.25. 


<AD  MEAN  CBD  STANDARD  DEVIATION 
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The  curves  presented  are  for  Kq  =  0.025,  0.10,  0.25  and  0.50  and  N  =  64, 
128,  and  512.  Values  of  the  mean  and  standard  deviation  for  other  values 
of  Kq  and  N  can  be  obtained  by  interpolating  from  the  curves 
presented. 


3.3  Statistics  of  the  Fourier  Coefficient  for  Non-Stationary 

Periodic  Input 

The  matched  filter  detector  was  described  in  Chapter  II  and  shown 
in  Figure  2-5.  The  first  function  performed  by  the  matched  filter  detector 
is  a  Discrete  Fourier  Transform  of  the  average  phase  estimate  4>a Cm) .  The 
mean  and  standard  deviation  of  the  average  phase  estimate  are  known.  The 
standard  deviation  is  a  function  of  the  signal  phase  which  changes  with 
time.  Therefore,  the  input  to  the  DFT  must  be  considered  to  be  a  non¬ 
stationary  process. 


The  real,  discrete  time  sequence  <j>  (m)  is  a  non-stationary  random 

a 


process  with  mean,  0  (m) ,  and  variance,  o.  (m) .  The  mean  and  variance  are 

a  <pa 

known  for  values  of  m  between  zero  and  (N_-l).  Over  this  period,  the  mean 

r 

and  variance  of  <J>a(m)  can  be  expanded  by  a  Discrete  Fourier  Transform, 


<j>  (m) 
a 


V1 

—  I 

nf  L 
k=0 


.  27Tkm 
C  (k)eJ  N  , 


m 


0<  m  <N  -1 

r 


(3-43) 


or 


V1 


4>  (m) 
a 


—  I 

N_  L 


\  2tt  km 

a  (k)  cos  — —  + 
m  N 

F 


k-0 


V1 


—  T 

nf  L 


i  /  \  \  i  2tt  km 

b  (k)  sin  - - 

m 


0£m<  N  p-l 


(3-44) 


k-0 
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where 


VI 


c  (k) 

in 


-  I 


.  2T!'km 


4>a(m)  e"J  Np  ,  0<k<  Nf-1 


m=0 


nf'1 


am(k>  =  Z  Vm)  COS  N^’  °-k-  V1 

F 

m=0 


Vi 


bm(k)  =  Z  Cm)  sin  0<k<NF-l 


m=0 


and 


Np-1 

2  i  .  2Trkm 

a<Jia  =  ?  Z  Cv(k)  eJ  Nf  ,  0<m<  Hp-1 

k=0 


>  (3-45) 


/ 


(3-46) 


or 


V1 


v(m) 


1  r*  >  2trkm 

•  Nl  Z  V10  cos  ~ 


K=0 


V1 


+  f  E  Vk)  sin  it5  •  °irai  v1 


k-0 


(3-47) 
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where 


C  (k) 

v 


2  2iTkm 

E  %(»>  e"3  » 


,  0<k<  Np-1 


a  (k) 
v 


Za  (m)  cos  N  ,  0<k<  N  -1 
<pB  - F 


(3-48) 


Vk>  =  E  °J(-0  sin  »  0£kl  *»_-! 


Sequential  samples  of  the  random  sequence  are  assumed  to  be  independent. 

The  discrete  time  sequence,  0  (m) ,  can  also  be  expanded  bv  a  Discrete 

a 

Fourier  Transform  over  the  period  Np.  It  is  necessary  to  know  the  moments 

of  the  resulting  Fourier  coefficient;  c.(k),  a,(k),  and  b,(k);  to  establish 

<P  9  <P 

the  statistical  performance  of  the  matched  filter  detector.  These  Fourier 
coefficients  are: 


CaOO 


.  2TTkm 

£  <b  (m)  e~3  N  ,  0<k<  N  -1 


Vk) 


£  0  (m)  cos  ■  ,  CKk£  N_-l 


(3-49) 


b.(k) 


Z.  ,  .  .  4TTkm  ^  , 

<P  (m)  sin  — —  ,  0<k<  N_-l 
a  Np  - F 


v 
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3.3.1  Mean  Value  of  the  Fourier  Coefficient 

The  mean  value  of  a  linear  summation  of  random  variables  is  equal 
to  the  linear  summation  of  the  mean  values  of  the  random  variables.  There¬ 
fore,  the  mean  values  of  the  Fourier  Coefficients  in  equation  (3-49)  are 
equal  to  the  right-hand  side  of  equation  (3-45).  Thus 

"  c»(k)  ) 

V0  ‘  a-<k)  i  <3'50> 

V«  ■  Vk> 

The  mean  of  the  Fourier  coefficients  for  the  DFT  representation  of  5  (n) 
is  equal  to  the  respective  Fourier  coefficients  of  the  DFT  expansion  of  the 
mean  4>  (m)  of  the  random  sequence. 

3.3.2  Variance  of  the  Fourier  Coefficients 

The  variance  of  a  linear  summation  of  uncorrelated  random  variables 
is  equal  to  the  summation  of  the  squares  of  the  linear  coefficients  times 
the  variance  of  each  random  variable.  Therefore,  variance  of  the  Fourier 
coefficients  is  given  by: 
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Comparing  the  above  equations  with  Equations  (3-48)  ,  the  covariance 
expressions  become 


Uc(k,£)  =  cv(k-A) 


y  (k,£)  =  |  a  (k-£)  +  |  a  (k+£) 

a  2  v  2  v 

y,  (k,£)  =  \  a  (k-A)  -  |  a  (k+A) 

D  Z  V  Z  V 


(3-59) 


y  (k,A)  =  ^  b  (k+A)  -  |  b  (k-A) 

ab  2  v  2  v 


3.3.5  Distribution  of  the  Fourier  Coefficients 

The  sine  and  cosine  Fourier  coefficients  are  formed  by  a  linear 

transformation  of  random  variables.  Therefore,  the  Fourier  coefficients 

are  also  random  variables  and  they  have  a  Gaussian  limiting  distribution. 

For  the  cases  under  consideration  it  will  be  assumed  that  N  is  large 

enough  that  the  central  limit  theorem  applies  and  the  distribution  of 

the  sine  and  cosine  Fourier  Coefficients  can  be  considered  to  be  Gaussian 

with  mean  and  variance  as  previously  derived. 

The  moments  of  the  Fourier  coefficients  for  the  real  discrete, 

random  sequence  $a(m)  are  summarized  in  Table  3-1.  If  the  Fourier 

Coefficients  of  the  DFT  expansion  of  the  variance  of  sequence  <j>  (m)  has 

a 

a  DC  term  which  is  large  with  respect  to  any  of  the  higher  order  coeffi¬ 
cients,  the  input  sequence  can  be  considered  to  be  stationary.  For  a 
stationary  sequence,  the  mcnents  given  in  Table  3-1  can  be  simplified  as 


TABLE  3-1 
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given  in  Table  3-2. 


Expanding,  the  variance  of  becomes 
K 

°x2  -  £  [ag2(k)  °a2<k)  +  bs2(k>  ab2W  +  2as<k)  Vk)  ^ab(k)] 

k-1 

K  K 

+  Z  Z  ta3(k)  as(Z,)ya(k’Jl)  +  bs(k)  bs(J')  Vk,5,) 

k-1  1-1 

l+k 

+  2  ag(k)bg(fc)  ab<k,fc)] 


(3-61) 


TABLE  3-2 


APPROXIMATE  VALUE  FOR  MOMENTS  OF  FOURIER  COEFFICIENTS 
ASSUMING  STATIONARY  INPUT  SEQUENCE 
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where  the  values  of  the  variance  and  covariance  terms  can  be  found  in 
Tables  3-1  or  3-2.  Since  the  Fourier  coefficients  are  assumed  to  be 
Gaussian,  the  input  to  the  threshold  detector  is  also  Gaussian. 


1  r  <XT  - 

P<V  =  7=“  exp  "  2'  I 

STn  L  2 Oyj  J 

The  probability  that  will  exceed  the  threshold  is 


(3-62) 


P(W 


-ifi-e  t  fl^l] 

L  '^x  'J 


(3-63) 


The  likelihood  functions  can  be  obtained  from  the  probability  desnity 
function  of  A„ 


poaT)  -  pttT  i  at  <  xT  +  dxT  |ho) 


(3-64) 


pl(xt)  =  P(XT  <  T  <  xT  +  dXT  |h1) 


(3-65) 


3.5  Receiver  Operating  Characteristics 

The  receiver  operating  characteristics,  ROC,  curves  relate  the 
probability  of  detection  versus  the  probability  of  false  alarm  for  a  given 
signal-to-noise  ratio.  These  curves  can  be  determined  from  the  properties 
of  the  desired  signal,  the  signal  phase,  the  properties  of  the  correlated 
noise  and  the  signal-to-noise  ratio.  The  likelihood  ratios  can  be  found 
in  the  following  manner: 

1.  Determine  the  input  noise  correlation  coefficients  p„  and  A  : 

N  N 

the  signal-to-noise  ratio  h  and  the  phase  <£g(t),  of  the 
desired  signal. 
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2.  Calculate  the  observable  signal  correlation  coefficient 
magnitude  Kq  and  phase  4>0  using  Equations  (3-12)  and  (3-13). 
Values  of  Kq  and  <J>Q  should  be  calculated  over  the  total  processing 
interval  (or  sequence  period) . 

3.  Determine  the  mean  and  variance  of  the  average  phase  over  the 
total  processing  interval.  NOTE:  it  is  assumed  the  averaging 
time  is  short  enough  that  the  signal  phase  can  be  considered 
constant  during  the  averaging  interval. 

4.  Calculate  the  Fourier  coefficients  for  the  mean  and  variance  of 
the  average  phase. 

5.  Calculate  the  statistical  moments  of  the  DFT  coefficients  using 
Table  (3-1)  [Table  (3-2)  if  the  average  phase  variance  is  inde¬ 
pendent  of  time] . 

6.  Determine  the  mean  and  variance  of  the  threshold  detector  input 
using  Equations  (3-60)  and  (3-61). 

7.  The  likelihood  function  is  calculated  using  Equation  (3-62). 

For  the  Hq  hypothesis,  the  mean  and  variance  of  the  average  phase 
will  be  constant  over  the  processing  period.  The  input  to  the  matched 
filter  is  stationary  and  the  values  for  the  moments  of  DFT  coefficients 
are: 


c(k) 


,  k  *  0 


_  2„  „  _  2„  %  NF  °ia 

aa  (k)  "  ab  00  =  ~ 2 


(3-66) 


(3-67) 


where 
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0^a  -  variance  of  average  phase  estimate  given  Hq 

Np  -  number  of  points  in  DFT. 

The  mean  and  variance  of  the  threshold  detector  input  are: 


K  =0 

To 


V  =  W  £  [a/(k)  +  bsZ(k)] 

k=l 


and  the  likelihood  ratio  is: 


P°(V  ‘  ^  “p 


2aXoZ  ' 


The  probability  of  false  alarm  is: 


fa  2 


■ 

1  -  e  r  f 

H 

Nil 

Q 

o 

(3-68) 


(3-69) 


(3-70) 


(3-71) 


The  likelihood  ratio  for  the  hypothesis  can  be  found  using 
the  steps  given.  The  probability  of  detection  can  then  be  determined 
from  the  p^(A^)  likelihood  function. 


57 


CHAPTER  IV 

DETECTION  OF  KNOWN  SIGNALS 

The  preceding  chapters  described  a  receiver  to  detect  a  dual 
channel  differential  phase  modulated  signal  and  developed  the  theoretical 
results  necessary  to  determine  the  receiver  operating  characteristic 
curves.  In  this  chapter  we  will  apply  those  results  to  the  detection  of 
"known  signals."  "Known  signals"  is  used  to  mean  signals  as  described  in 
Chapter  II  with  a  known  signal  phase  delay  function.  One  example  will 
be  examined.  The  example  is  a  linear  delay  which  is  often  encountered  in 
a  practical  application.  Radiation  from  a  source  moving  with  a  constant 
velocity  will  result  in  an  approximate  linear  phase  delay  between  the 
output  of  two  spatially  separated  sensors.  The  matched  filter  configura¬ 
tion  is  determined,  the  mean  and  variance  of  the  average  phase  estimate 
calculated  in  both  the  time  domain  and  the  frequency  donian,  and  the  ROC 
curves  developed.  The  average  phase  estimate  is  calculated  for  both 
the  angle  average  and  vector  average  cases  so  that  the  effect  of  the 
averaging  method  on  the  receiver  operating  characteristics  can  be  examined. 

Detection  of  these  known  signals  is  evaluated  under  varying  noise 
correlation  conditions.  The  magnitude,  K^,  and  phase,  ,  of  the 
noise  correlation  coefficients  considered  are: 
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Two  cases  are  investigated  using  the  preceding  correlated  noise  conditions. 
For  the  first  case,  the  noise  correlation  is  known  a  priori  allowing  the 
detection  threshold  to  be  fixed  for  a  desired  false  alarm  rate.  In  the 
second  case,  the  noise  correlation  is  unknown  a  priori  and  the  threshold 
is  set  based  on  an  assumed  correlation  condition.  The  analytical  ROC 
curves  are  generated  for  two  threshold  settings. 

The  examples  considered  have  6400  points  in  the  processing  interval 
and  each  interval  represents  one  period  of  the  time  delay  function.  The 
6400  points  are  processed  by  averaging  100  points  and  performing  a  64-point 
DFT.  No  attempt  is  made  to  investigate  an  optimum  trade-off  between  the 
number  of  points  averaged  and  number  of  points  in  the  DFT. 


4,1  Linear  Delay  With  Known  Noise  Correlation 


The  linear  (or  ramp)  phase  delay  function  is: 


4>  (t)  =  d>  (-1  +  ?t_) ,  0  <  t  <  T  . 

s  P  T  ~ 
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The  phase  delay  of  Che  sampled  signal  becomes: 

$  (n)  =  ij>  (-1  +  2H) ,  0  <  n  <  N-l  .  (4-1) 

s  P  N  ~  ~ 

For  the  example  being  investigated  <J>  was  chosen  to  be  tt/2  radians 

P 

(90°). 

The  Fourier  transform  of  the  linear  delay  is: 


(4-2) 

(4-3) 


The  configuration  used  in  implementing  the  filter  is  an  approximation 
of  the  ideal  matched  filter  in  that  only  the  first  three  harmonics  are 
used  as  shown  in  Figure  4-1.  The  configuration  was  chosen  to  simplify 
implementation  and  because  after  the  averaging  process  the  higher  harmonics 
fall  off  faster  than  the  inverse  of  k. 

The  likelihood  functions  are  again  found  using  the  procedure  from 
Chapter  III.  The  steps  will  not  be  presented  in  detail.  Any  approxima¬ 
tions  will  be  discussed  and  the  results  presented. 

Figures  4-2  through  4-4  show  the  correlation  coefficient  between 
the  observable  signals,  the  angle  average  mean  and  variance  and  the  vector 
average  mean  and  variance,  respectively.  These  variables  are  all  shown  for 
a  signal-to-noise  ratio  of  -12  dB . 

The  mean  and  variance  of  the  input  to  the  threshold  detector  are 
based  on  the  first  three  harmonics.  The  mean  of  A  is: 


FIGURE  4-2. 


Variation  in  K  and  i  for  Linear  Delav  and  SNR 
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=  ±r  [-b  (1)  -  \  b  (2)  -  ^  b  (3)] 
64  m  2  m  3  m 


(4-4) 


and  the  variance  can  be  calculated  assuming  the  average  phase  estimate 
is  stationary  with  a  variance  equal  to  the  average  of  the  variance  of 
the  individual  estimates  over  a  processing  interval. 


a 


\ 


1 

2(64) 2 


~  c  (0) 
36  v 


(4-5) 


Equation  (4-5)  was  used  in  calculating  the  ROC  curves.  In  Figures 

_3 

4-5  and  4-6,  the  ROC  curves  for  a  probability  of  false  alarm  of  10 
and  all  noise  correlation  conditions  are  combined  for  angle  average  and 
vector  average  phase  estimates,  respectively.  A  comparison  of  the 
receiver  performance  as  a  function  of  noise  correlation  can  be  made.  The 
variation  in  signal-to-noise  ratio  for  a  given  p^  is  less  than  1.2  dB 
for  the  angle  average  case.  The  vector  average  method  results  in  receiver 
performance  which  is  more  sensitive  to  noise  correlation. 

The  ROC  curves  for  the  angle  average  phase  estimate  case  can  be 
further  understood  with  the  aid  of  Figure  4-3.  The  waveform  shape  of 
the  mean  of  the  average  phase  estimate  remains  highly  linear  over  all  the 
noise  correlation  conditions.  The  variance  varies  over  a  3:1  range  as 
the  noise  correlation  changes  but,  in  general,  the  peak-to-peak  value  of 
the  mean  increases  as  the  variance  increases.  Thus,  the  detection  per¬ 
formance  when  using  angle  average  phase  estimate  is  relatively  insensitive 
to  the  noise  cor 'elation. 


The  detection  performance  when  using  the  vector  average  phase 
estimate  can  be  explained  with  the  help  of  Figure  4-4.  Consider  first, 


FIGURE  A- 5.  Detection  Performance  for  Linear  Delay  and  Angle  Average  Phase  Estimate 
Noise  Correlation  Known  a  Priori. 


FIGURE  4-6.  Detection  Performance  for  Linear  Delay  and  Vector  Correlation  Known  a  Priori 
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the  case  of  correlated  noise,  =  0.2,  =  it.  The  variance  of  the 

phase  estimate  is  2.86  when  no  signal  is  present.  The  standard  deviation 

of  the  test  statistic  under  the  Hq  hypothesis  is  0.295.  At  the  high 

signal-to-noise  ratio,  the  maximum  peak  value  of  the  mean  phase  estimate  is 

ir/2.  Thus,  the  maximum  mean  value  of  the  test  statistic  is  0.458.  In 

-3 

order  to  obtain  a  probability  of  false  alarm  of  10  ,  the  threshold 

must  be  set  3.09  times  the  standard  deviation  or  greater  than  maximum 

mean  value  of  the  test  statistic.  As  the  signal-to-noise  ratio  increases 

the  variance  and  mean  decreases  and  the  probability  of  the  test  statistic 

exceeding  the  threshold  decreases  at  signal-to-noise  ratios  greater  than 

_2 

-6  dB.  For  probability  of  false  alarm  greater  than  6  X  10  ,  the  maximum 

mean  value  of  the  test  statistic  will  exceed  the  threshold.  This  behavior 
could  be  avoided  by  increasing  the  number  of  points  in  the  matched  filter 
DFT. 

Examination  of  mean  phase  estimate  in  Figure  4-4  shows  the  waveform 

is  badly  distorted  for  noise  correlation  of  K  =  0.5,  <J>  =  tt/2.  With 

N  N 

noise  correlation  phase  of  tt/2  and  signal-to-noise  ratio  less  than  the 
noise  correlation  magnitude,  the  mean  vector  average  phase  estimate 
waveform  is  no  longer  matched  to  the  filter  in  the  detector.  As  the 
signal-to-noise  increases  the  mean  vector  average  phase  estimate  waveform 
ap 'in  becomes  approximately  linear.  This  noise  ratio  decreases  in  Figure 
4—6  when  the  noise  correlation  phase  is  tt/2. 


I 


4. 2  Unknown  Noise  Correlation,  Fixed  Threshold  Detector 

Previously,  it  was  assumed  that  the  magnitude  and  phase  of  the 
noise  correlation  coefficient  were  known  and  the  detector  threshold  could 
be  set  for  a  given  probability  of  false  alarm.  It  is  also  of  interest  to 
examine  the  case  of  a  fixed  threshold  when  the  noise  correlation  is  un¬ 
known  . 

Figures  4-7  and  4-8  are  the  receiver  operating  characteristics  when 

the  angle  average  method  is  used  to  estimate  the  average  phase.  The  thres- 

-3 

hold  was  set  to  give  a  probability  of  false  alarm  of  10  ;  in  Figure  4-7 

the  threshold  was  set  based  on  uncorrelated  noise  and  in  Figure  4-8  a 
noise  correlation  coefficient  of  =  0.5,  0^=0  was  assumed.  These 
curves  can  be  compared  directly  with  Figure  4-5  where  the  noise  correlation 
was  assumed  known. 

For  the  values  of  actual  noise  correlation  considered  in  Figure  4-7, 

_3 

the  worse  case  probability  of  false  alarm  would  increase  to  2.4  x  10  at 

a  noise  correlation  of  =  0.2,  =  tt.  For  an  actual  noise  correlation 

of  =  0.5,  =  0  or  =  0.7,  <|)  ~  ir/ 4  the  receiver  sensitivity  suffers 

a  1.1  dB  loss.  In  Figure  4-8,  the  receiver  sensitivity  suffers  no  loss, 

_2 

but  the  probability  of  false  alarm  increases  to  1.9  x  10  if  the  actual 

noise  correlation  were  =  0.2,  4>N  =  tt,  the  probability  of  false  alarm 
-3 

would  not  exceed  10  for  any  of  the  actual  cases  considered  and  the  loss 
of  receiver  sensitivity  would  be  1.8  dB.  Thus,  the  performance  of  the 
receiver  is  relatively  insensitive  to  noise  correlation  coefficient  when 
using  an  angle  average  phase  estimate  and  for  the  values  of  correlation 


considered. 


Detection  Performance  for  Linear  Delay,  Angle  Average  Phase  Estimate 
Unknown  Noise  Correlation,  and  Threshold  Set  Assuming  K  =  0,5,  <f>  = 
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The  receiver  operating  characteristics  when  the  vector  average 
phase  estimate  is  used  are  shown  in  Figures  4-9  and  4-10  for  the  same 
cases  as  considered  previously.  In  Figures  4-9  and  4-10  the  threshold 
was  set  assuming  noise  correlation  of  =  0.0  and  =  0.5,  =  0, 

respectively.  These  receiver  operating  characteristics  are  sensitive  to 
actual  noise  correlation  coefficient.  In  Figure  4-9,  the  probability  of 

_3 

false  alarm  increases  to  2.5  x  10  for  an  actual  noise  correlation  of 
=  0.2,  0  =  it.  For  other  noise  correlation  conditions,  the  receiver 

sensitivity  decreases  such  that  for  a  signal-to-noise  ratio  of  -4  dB  the 

-4 

probability  of  detection  is  less  than  10  .  The  curves  in  Figure  4-10 

show  a  worse  case  decrease  in  expected  receiver  sensitivity  of  6.8  dB 
while  the  probability  of  false  alarm  increases  to  0.45  when  the  noise 
correlation  is  =  0.2,  c|)^  =  tt.  It  would  be  undesirable  to  use  the 
vector  average  phase  estimate  if  the  noise  correlation  coefficient  were 
not  exactly  known. 

4 . 3  Comparison  with  Cross  Correlation  Detector 

The  cross  correlation  receiver  is  commonly  used  to  detect  Gaussian 
signals  in  Gaussian  noise.  To  detect  a  dual  channel,  differential  phase 
modulated  signal  in  correlated  noise,  the  cross  correlator  receiver  in¬ 
cludes  a  noise  decorrelator  or  spatial  pre-whitener  and  phase  compensation. 
The  detection  performance  of  this  cross  correlation  receiver  is  compared 
with  the  phase  matched  filter  receiver  for  detection  sensitivity  and  the 
sensitivity  of  the  receivers  to  the  noise  correlation  coefficient. 


The  cross  correlation  detector  with  spatial  prewhitening  and  signal 
phase  compensation  is  described  in  Appendix  B.  The  necessary  equations 


Noise  Correlation,  and  Thres 
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to  calculate  the  receiver  operating  characteristic  curves  are  also 

contained  in  Appendix  B.  The  ROC  curves  for  the  linear  time  delay  and 

-3 

a  probability  of  false  alarm  of  10  were  generated  for  two  cases: 

1)  the  noise  correlation  known,  and  2)  noise  correlation  unknown  but 
threshold,  prewhitening  and  phase  compensation  set  assuming  a  noise  corre¬ 
lation  value.  These  ROC  curves  are  compared  with  corresponding  curves 
for  the  signal  phase  matched  filter  detector  using  either  angle  average 
or  vector  average  phase  estimation. 

4.3.1  Known  Noise  Correlation 

The  ROC  curves  for  a  cross  correlation  detector  using  spatial 
prewhitening  and  signal  phase  compensation  are  shown  in  Figure  4-11.  These 
curves  represent  the  case  when  the  threshold  was  set  to  give  a  probability 

-3 

of  false  alarm  of  10  and  thus  correspond  to  the  ROC  curves  in  Figures 
4-5  and  4-6  for  the  receiver  using  a  signal  phase  matched  filter.  For 
the  values  of  noise  correlation  considered,  '  <e  cross  correlation  receiver 
using  spatial  prewhitening  is  3.3  to  4.7  dE  more  sensitive  than  a  receiver 
using  angle  average  phase  estimate  and  a  phase  matched  filter  detector. 

The  cross  correlation  receiver  is  0.7  to  8.8  dB  more  sensitive  than  a 
receiver  using  vector  average  phase  estimate  and  a  phase  matched  filter 
detector. 

4.3.2  Unknown  Noise  Correlation 

Two  cases  were  examined  for  the  cross  correlation  receiver  with 
the  spatial  prewhitening,  signal  phase  compensation  and  threshold  set 
assuming  a  given  noise  correlation  and  the  actual  noise  correlation 
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being  a  different  value.  The  ROC  curves  in  Figures  4-12  and  4-13  are  for 
assumed  noise  correlation  of  =  0.0  and  =  0.5,  <t>^  ■  0,  respectively. 

If  the  assumed  noise  is  uncorrelated  and  the  actual  noise  correlation  is 
Kjj  -  0.5,  <j>N  =  0;  *  0.5,  <)>N  =  ir/2  or  ■  0.7,  <J>N  =*  it/4,  the  probability 

of  detection  will  be  greater  than  0.9999.  However,  the  probability  of 
false  alarm  will  also  be  greater  than  0.9999.  The  same  applies  if  the 
assumed  noise  correlation  is  =  0.5,  4>N  *  0  and  the  actual  noise 
correlation  is  =  0.2,  4>n  *  tt.  To  reduce  the  probability  of  false  alarm, 
it  would  be  necessary  to  increase  threshold  and  suffer  the  corresponding 
loss  in  sensitivity. 

The  ROC  curves  in  Figures  4-12  and  4-13  can  be  compared  directly 
with  Figures  4-7  through  4-10  for  the  receivers  using  signal  phase  matched 
filter  detector.  The  receiver  which  is  least  sensitive  to  the  actual 
noise  correlation  is  the  receiver  using  angle  average  phase  estimation 
and  a  signal  phase  matched  filter. 


and  Threshold  Set  Assuming 
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CHAPTER  V 

SUMMARY  AND  CONCLUSIONS 


5.1  Summary 

This  research  studied  the  signal  detection  characteristics  of  a 
receiver  with  two  input  channels.  The  signals  received  in  both  channels 
are  random,  narrowband  processes  which  are  identical  except  for  a  known, 
variable  channel-to-channel  time  delay.  The  known  time  delay  is  represented 
by  a  carrier  frequency  phase  shift  when  expressing  the  narrowband  random 
signals.  The  noise  is  also  a  narrowband,  random  process.  There  can  be 
correlation  between  the  noise  in  each  of  the  input  channels. 

Chapter  II  explained  the  generalized  received  signal  and  noise 
model  and  described  the  receiver  being  studied.  The  receiver  estimates 
the  average  phase  difference  between  the  observable  signals  in  each  of 
the  input  channels.  Methods  of  forming  the  average  phase  estimate  consid¬ 
ered  are  angle  averaging  and  vector  averaging.  The  average  phase  signal  is 
passed  through  a  filter  whose  transfer  function  is  matched  to  the  known 
time  varying,  signal  phase.  The  filter  is  implemented  using  a  discrete 
Fourier  transform.  The  output  of  the  filter  is  threshold  detected  to 
determine  the  presence  of  a  signal  with  the  known  phase  shift. 

The  theoretical  statistics  required  to  determine  the  detection  per¬ 
formance  of  the  receiver  were  developed  in  Chapter  III.  The  statistics 
are  functions  of  the  input  signal-to-noise  ratio,  the  correlation  of  the 
noise  in  the  two  receiver  input  channels,  the  period  of  the  phase  averages 
and  the  number  of  points  used  in  the  discrete  Fourier  transform.  The 
moments  of  the  vector  average  phase  estimate  were  derived.  The  variance  of 


the  average  phase  signal  varies  with  time;  therefore,  the  moments  of  the 
Fourier  coefficients  were  derived  assuming  a  non-stationary  periodic  in¬ 
put.  The  receiver  likelihood  functions  were  determined  so  the  probability 
of  false  alarm  and  probability  of  detection  could  be  calculated. 

In  Chapter  IV,  the  detection  performance  was  determined  for  an  example 
of  possible  signal  phase  shift  function.  The  phase  shift  function  consider¬ 
ed  was  a  linear  phase  shift.  The  receiver  operating  characteristic  curves 
were  generated  for  several  noise  correlation  conditions  assuming  the  noise 
correlation  was  known  a  priori.  Finally,  receiver  operating  characteris¬ 
tic  curves  were  generated  when  the  detection  threshold  was  set  for  an 
assumed  noise  correlation  but  the  actual  noise  correlation  was  not  known 
a  priori. 

A  performance  comparison  of  the  receiver  using  an  average  phase 
estimate  and  phase  matched  filter  detector  to  a  cross  correlation  detector 
was  made.  The  cross  correlation  detector  included  a  noise  decorrelator 
and  phase  shift  compensator.  Comparisons  were  made  assuming  the  noise 
correlation  was  known  a  priori  and  then  assuming  the  threshold  was  set 
for  a  given  noise  correlation  but  the  noise  was  not  known  a  priori. 

5. 2  Conclusions 

It  was  stated  that  the  signal  processing  was  intended  to  detect 
at  low  signal- to-noise  ratios  and  over  a  variety  of  noise  correlation 
values.  The  receiver  proved  to  be  robust  under  these  conditions  when  the 
angle  averaging  method  of  estimating  the  average  phase  was  used.  For  the 
noise  correlation  conditions  investigated,  the  variation  in  signal- to-noise 
ratio  for  a  given  probability  of  detection  was  less  1.2  dB.  The 
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performance  when  using  the  vector  average  phase  estimate  was  more  sensitive 
to  the  value  of  the  noise  correlation.  When  the  noise  correlation  is  known 
exactly  a  priori,  the  choice  of  phase  averaging  method  depends  on  the 
exact  value  of  the  correlation  and  the  desired  probability  of  false  alarm. 

In  the  case  of  a  fixed  threshold  based  on  an  assumed  value  of  noise  cor¬ 
relation  where  the  actual  noise  correlation  is  not  known  a  priori,  the 
performance  using  the  vector  average  phase  estimate  is  unsatisfactory. 

The  variance  of  the  test  statistics  and  hence  the  threshold  setting  is 
sensitive  to  the  actual  value  of  the  noise  correlation.  The  angle  average 
phase  estimate  is  less  sensitive  to  the  actual  value  of  the  noise 
correlation.  When  the  noise  correlation  is  not  known  a  priori  and  the 
angle  average  phase  estimate  is  used,  a  threshold  setting  can  be  found 
which  will  give  variation  in  signal-to-noise  ratio  of  less  than  6  dB  for 
a  given  probability  of  detection  and  a  probability  of  false  alarm  below  a 
desired  level.  For  the  noise  correlation  values  studied,  the  variation 
in  signal-to-noise  ratio  for  50%  probability  of  detection  was  2  dB. 

A  comparison  of  a  cross  correlation  receiver  with  noise  decorrelator 
and  phase  compensator  and  a  receiver  using  an  average  phase  estimate  and 
phase  matched  filter  detector  was  made  for  the  noise  correlation  known  a 
priori  and  for  a  threshold  based  on  an  assumed  noise  correlation  when  the 
actual  correlation  was  not  known  a  priori.  The  cross  correlation  receiver 
is  approximately  4  dB  more  sensitive  than  a  receiver  using  an  angle 
average  phase  estimator  when  the  noise  correlation  is  known  a  priori.  The 
cross  correlation  receiver  is  1  to  9  dB  more  sensitive  than  a  receiver 
using  an  vector  average  phase  estimator  when  the  noise  correlation  is  known 
a  priori.  If  the  noise  correlation  is  not  known  a  priori,  the  performance 
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of  the  cross  correlation  detector  is  unsatisfactory. 

The  major  accomplishment  of  this  research  is  the  development  of  a 
receiver  with  fixed  detection  threshold  and  signal  processing  which  will 
detect  a  differential  phase  modulated  signal  in  correlated  noise  when  the 
value  of  the  noise  correlation  is  unknown.  The  phase  matched  filter 
receiver  with  an  angle  average  phase  estimator  can  be  used  over  a  wide 
range  of  noise  correlation  values.  Although  it  is  not  as  sensitive  as  the 
cross  correlation  receiver  under  known  noise  conditions,  its  sensitivity 
does  not  degrade  as  much  when  the  noise  correlation  is  not  known  a  priori. 

5.3  Future  Research 

Several  areas  need  further  research  to  fully  understand  and  optimize 
a  receiver  to  detect  a  dual  channel,  differential  phase  shifted  signal. 
Some  of  these  topics  are  discussed  briefly. 

The  detection  of  known  signals  was  performed  with  a  given  phase 
averaging  period  and  a  given  number  of  points  in  the  discrete  Fourier 
transform.  Mo  attempt  was  made  to  optimize  these  variables  or  the  trade¬ 
off  between  them.  A  trade-off  study  should  be  performed  to  determine  the 
total  processing  period  and  the  division  of  the  samples  within  the  process¬ 
ing  period  between  the  phase  averages  and  discrete  Fourier  transform. 

The  detection  of  known  signals  assumed  a  peak  phase  for  the  desired 
signal.  In  some  systems,  the  magnitude  of  the  peak  phase  may  be  a  design 
parameter.  The  curves  of  the  mean  average  phase  have  a  peak  at  a  value  of 
signal  pnaae  that  is  a  function  of  the  signal-to-noise  ratio  and  noise 
correlation  coefficient.  The  presence  of  this  peak  implies  that  there  may 
be  a  value  peak  signal  phase  that  is  optimum. 
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APPENDIX  A 

SELECTED  MIXED  CENTRAL  MOMENTS  OF  QUADRI VARIATE 
GAUSSIAN  DENSITY  FUNCTION 


The  mean,  variance,  and  covariance  of  the  real  and  imaginary 
components  of  a  product  of  two  complex  variables  can  be  found  from  the 
appropriate  mixed  central  moments  of  the  joint  density  function  of  the 
components  of  the  two  complex  variables.  The  product  of  a  complex  variable 
multiplied  by  the  complex  conjugate  of  the  second  complex  variable  is 
given  by 


uv*  ■  (u  +  ju  )  (v  -  jv  ) 
c  J  s  c  J  s 


uv*  =*  (uv  +  uv)  +  j(uv  -  uv), 
c  c  s  s  sc  cs 


(A-l) 


The  real  and  imaginary  components  are: 


x  =  Re  Tuv*l  =  u  v  +  u  v 
J  c  c  s  s 


(A-2) 


y  3  Im  [uv*l  =  u  v  -  u  v 
J  l  J  s  C  C  X 


(A-3) 


The  mean,  variance,  and  covariance  values  can  be  obtained  from  the 

mixed  central  moments  of  the  joint  probability  distribution  of  the 

variables  u  ,  u  ,  v  ,  and  v  • 
esc  s 

The  variables  u  ,  u  ,  v  ,  and  v  are  zero  mean  Gaussian  variables, 
c  s  c  s 

described  by  covariance  matrix: 


R  =*  [E  (Z  Z'}] 


(A-4) 


where  the  vector  Z  is  related  to  the  variables  uc ,  ug ,  vc,  and  vg  by 
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(A-5) 


The  ij  th  term  of  R  is  equal  to  the  jith  term  so  that  the  covariance 
matrix  is  symmetric.  The  quadrivariate  Gaussian  density  function  is 
given  by 

p(i>’  HrE'li 


and  the  characteristic  function  is 


C(j^) 


=  e 


Rw 


or 


4 

C(jw)  -  exp  {-  y  Y. 

i-1 


4 

Z  “k  rki  V 

k=l 


(A-7) 


(A-8) 


where 


r 


ij 


id- 


E  {ZjZ  } 

w j »  cj  ^  ] 


bi  b2  b3  bi» 

The  mixed  moment  E{x^  X2  x^  x^  }  where  b^  are  positive  integers 


can  be  found  from 
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(A-ll) 
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2 

Similarly  for  Eiz^^  z^  z^},  the  only  non-zero  term  are  of  the  form 


u) .  iii .  0), 
i  3  k 


.  therefore, 


E{w  w.  w  }  =  r..  r..  +  2  r .  .  r., 

i  J  k  ii  jk  ij  ik 


(A-12 ) 


and 


E^Z1Z2Z3Z4^  *  r12r34  +  r13r24  +  r14r23 


(A-13) 


The  terms  of  interest  are: 


E{u  v  } 
c  c 

= 

r13 

E<Vs} 

= 

r14 

E{usvc} 

a 

r23 

E{u  v  } 
s  s 

= 

r24 

E{u  2v  2 
c  c 

} 

=  rll 

r33 

+ 

2 

2 

r13 

E(u  2v  2 
c  s 

} 

rll 

r44 

+ 

2 

2 

r14 

2  2 
E{u  v 
s  c 

} 

r22 

r33 

+ 

2 

2 

r23 

_  r  2  2 

Elu  v 
s  s 

} 

r22 

r44 

+ 

2 

2 

r24 

E{u  2v  v 
c  c 

s} 

=  rll 

r34 

+ 

2 

r13  r14 

2 

E{u  v  v 
s  c 

S} 

r22 

r34 

+ 

2 

r23  r24 

E{u  u  v 

CSC 

2) 

r33 

r12 

+ 

2 

r13  r23 

E{u  u  v 
CSS 

2} 

r44 

r12 

+ 

2 

r14  r24 

(A-14) 


Using  equations  (A-13)  and  (A-14) ,  the  moments  of  the  real  and  imaginary 
components  are: 


r,,  +  r 
13  24 


(A-15) 
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APPENDIX  B 

NOISE  DECORRELATOR,  SIGNAL  PHASE  COMPENSATION, 

CROSS  CORRELATION  RECEIVER 

The  cross  correlation  receiver  is  shown  in  Figure  B-l.  The 
observable  signals  are  processed  by  a  noise  decorrelator  and  phase 
compensator.  The  signal  processing  is  then  split  into  two  channels. 

The  first  channel  estimates  the  covariance  and  the  second  channel 
estimates  the  power  level.  The  weighted  differences  of  the  two  process¬ 
ing  channels  is  formed. 

If  the  weighted  difference  is  positive,  the  signal  present  hypothesis 
is  selected. 


Noise  Decorrelator 

The  transfer  function  of  the  noise  decorrelator  is: 


1 


0 


“d 


-S  _ L_ 

/l'h2 


(B-l) 


where  and  <j>^  are  the  estimated  magnitude  and  phase  of  the  noise 
correlator  and  coefficient.  The  covariance  matrix  of  the  input  signals 
is : 


1 


R  - 

r 


K 


o 


J*o 

e 
o 


( B—  2 ) 
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Phase  Compensation 

The  phase  compensation  operates  on  the  observable  signal  to  cause 
the  covariance  terms  to  be  real.  This  is  accomplished  with  the  transfer 
function: 
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The  signal  phase  is  a  time  varying  function.  The  time  delay  compensation 

then  must  also  be  a  function  of  time.  If  R  is  put  in  terms  of  the  real 

u 

and  imaginary  components  of  u  the  isomorphic  counterpart  for  the  general 
covariance  matrix  becomes: 


R  = 
u 


1 

0 


0 

1 


u 


0 


where 


1  +  s2  - 2  K<,s  “*(vv 


p  “Re 
u 


\  *  Im 

u 


J  j  <f»  -  j  <j> 

o  „  _  n*  _  ml 


<y  -v  > « 


^o*  3V  *j*»l 


(Koe  “Kjje  ")  e 


/  1-K 


-s 


(B-10) 


(B-ll) 


(B-12) 


(B-13) 


Covariance  Estimator 

The  covariance  estimator  performs  the  operation: 

N-l 

\  |  £  x(n)  (B-14) 

n-0 


where 
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x(n)  =  Re  (u^(n)  u2*(n) } 


The  mean  and  variance  of  x  can  be  found  from  the  expression  developed 
in  Appendix  A. 


x  =  2  a  p 
o  Ku 


2  o  2  2  ,  2s 

a  =  2  a  (a  +p  -  X  ) 
X  o  u  Ku  u 


The  x  variable  is  averaged  to  form  an  estimate  of  the  covariance. 
The  covariance  estimate  X^  is  Gaussian  distributed  with  a  mean  and 


(B-15) 


(B-16) 


(B-17) 


variance 


2  N-l 


X1 =  Hr  \  pu(n) 

n=*0 


(B-18) 


4  N-l 


J-  E  [au^n>  +  pu2<n)  *  Xu2^n)  ] 


(B-19) 


Power  Estimator 


The  power  level  estimate  is: 


X2  =  I  Z  ul(n)  ui*(n)  I  a2(n)  u2*(n) 


(B-20) 


The  power  level  estimator  is  assumed  to  be  Guassian  with  mean  and  variance: 


X,  *  2  o 
2  o 


1  +  i  E  au2(n) 
n-0 


(B-21) 
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N-l  \ 

I  [V<»>  +  2  »u2(n)  +  2  V2^])  • 

n=0  I 


(B-22) 


The  covariance  between  the  power  level  estimate  and  the  covariance 
estimate  is: 


8o  4  N_1 

WX  N~  J  I  [pu(n)  1  +  au2(n)] 

n«0 


(B-23) 


Test  Statistic 

The  test  statistic  is  the  weighted  difference  of  the  covariance 
estimate  and  the  power  level  estimate 


X  -  A,  -  G\_  . 

cl  2 


(B-24) 


The  mean  variance  of  the  test  statistic  is: 


N-l  N-l 

f  I  °u(">  -  G[1  +  l  I  “u2^] 


n“0 


n»0 


(B-25) 


4  N-l 


=  \{l  I  +pu2^>  -  \>>] 


n*0 


N-l 


+  G2  +  J-  E  [au2(n)  +  2  pu2(n)  +  2 


n“l 


N-l 


-  IT  £  [  1  +  «u2(n>] 


n«0 


(B-26) 


For  the  case  of  noise  only  and  the  noise  and  decorrelator  are 


perfectly  matched,  the  mean  and  variance  is: 


X  =  -  4o  ‘  G 
c  n  N 


8  a. 


N  1  +  G2 


The  signal  present  hypothesis  is  selected  if  the  test  statistic  is 
greater  than  zero.  The  probability  of  false  alarm  is  given  by 


fa 


erf 


-X  n 
c 


/2 


c  n 


fa 


*[■- 


erf 


I  /2N  G 


/r/z? 


The  probability  of  detection  is  given  by: 


fa 


»[- 


erf 


-X 


c  s+n 


J~1 


s+n 
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